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Abstract 

A canonical formalism of /(i?)-type gravity is proposed, resolving the prob- 
lem in the formalism of Buchbinder and Lyakhovich(BL). The new coordi- 
"q I nates corresponding to the time derivatives of the metric are taken to be 

its Lie derivatives which is the same as in BL. The momenta canonically 
conjugate to them and Hamiltonian density are defined similarly to the for- 
malism of Ostrogradski. It is shown that our method surely resolves the 
problem of BL. 

X' 
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1 Introduction 

Einstein gravity explains the observed universe fairly well. It has however some theoretical 
drawbacks, most important of which is the problem of initial singulariryfl, 2]. This problem 
is usually interpreted to imply the limit of applicability of Einstein gravity. There are 
two possible ways to approach this problem. One is to modify the theory of gravity in the 
classical framework. The most popular modified theory is the higher-curvature gravity(HCG) 
theories [3, 4, 5]. Another is to quantize gravity. Quantum gravity has a long history but 
its completion seems to require still a long time. 

HCG is also required in quantum field theory in curved spacetime[6], string perturbation 
theory [7]. Recently HCG is applied to cosmology to explain, e.g. inflation [8] or expansion of 
the present universe[10] which seems to be in the stage of an accelerated expansion[12, 13, 14]. 
In the former a term proportional to the square of the scalar curvature is added to the 
Lagrangian density of Einstein gravity[8, 9] and in the latter a term proportional to the 
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inverse of the scalar curvature is added[10]. In the first order Palatini formalism, a wider class 
of HCG's is also applied to cosmology to explain the present expansion of the universe. [11] 
Thus the reality of HCG is increasing. 

In this work we adress the problem in the canonical formalism of f(R) type HCG and 
propose a consistent formalism. The curvatures contain the second order time derivatives of 
the metric which in HCG cannot be removed by partial integration as in Einstein gravity. 
Thus HCG is a theory with higher order time derivatives(HDT). Standard procedure for 
the canonical formalism of HDT has been given by Ostrogradski[15]. In the usual method 
of canonical formalism of gravity in terms of ADM variables, this method however is not 
applicable directly. The reason is that the curvatures depend on the time derivatives of the 
lapse function and the shift vector so that they obey field equations leading to the breaking 
of general covariance. 

This problem is resolved by the method of Buchbinder and Lyakhovich(BL)[16, 17] which 
generalizes the choice of the new generalized coordinate corresponding to the time derivatives 
of the original generalized coordinates. Applied to the HCG, the new generalized coordinetes 
corresponding to the time derivatives of the metric was chosen to be the extrinsic curvature. 
The time derivatives of the lapse and shift are absorbed in the time derivatives of the extrinsic 
curvature. However, in this method, a change of the original generalized coordinates induces 
a change of Hamiltonian as will be explained below, which is not the case in theories without 
higer order time derivatives. For example, in the case of FRW spacetime, the Hamiltonian 
is different whether we use the scale factor or its logarithm as the generalized coordinate. 

We propose a canonical formalism which resolves this problem. We combine the advanta- 
geous points of both the method of Ostrogradski and BL. We choose the extrinsic curvature 
as the new generalized coordinate as in BL and define the momenta canonically conjugate to 
them similarly to Ostrogradski's method. The extrinsic curvature is the Lie derivative of the 
three metric which reduces to the time derivative in flat space. The resulting Hamiltonian 
is shown in fact invariant under the tarnsformations of the original generalized coordinates 
which do not change the 3-dimensional metric of the hypersurface of constant time. 

In section 2, we explain how the Hamiltonian changes under the transformation of the 
original generalized coordinate in the method of BL. Section 3 is devoted to the presentation 
of a new canonical formalism for a HCG of /(-R)-type. In section 4 we demonstrate the 
invariance of the Hamiltonian under the two kinds of transformations of the generalized 
coordinates, the three metric. Summary and discussion are given in section 5. 



2 A problem in the method of Buchbinder and 
Lyakhovich 

In this section we show that in the method of BL the Hamiltonian changes under the trans- 
fornation of the generalized coordinates using a simple model. Let us consider a system 
described by a Lagrangian 

L = L(q\q\<f) (2.1) 
in which case the new generalized coordinates are defined as 

Q l = q l . (2.2) 

In the method of BL, the Lagrangian is modified using the Lagrange multiplier method 
so that the definition (2.2) is derived from variational principle. Denoting the modified 
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Lagrangian as L*, it is given as 

L* = L + Pl (q l - Q l ), L — L(q\ Q\ Q*). (2.3) 

The multipliers Pi are the momenta canonically conjugate to q l . The canonical formalism 

is obtained by the Legendre transformation starting from L*. Momenta P l canonically 
conjugate to Qi are given by 

P< . 9 JC = a. (2.4) 

Then the Hamiltonian if* is given as 

#* = p.g* + p.Q* — L* — Vi Q l + PiQ* - L. (2.5) 

Now we make transformations of the generalized coordinates: 

V = fitf\ or q^g'tf). (2.6) 

Momenta conjugate to 0* are denoted as 7Tj, new generalized coordinates as $j and momenta 
canonically conjugate to them as Ilj. These variables are not uniquely related to old variables. 
Since Q % are related to q % through the variational principle, their relations to are not fixed 
a priori. The solutions are, however, Q % = q\ it is natural to assume the transformation of 
Q l to be the same as q\ Differentiation leads to 

f = £^ = Mji, (2.7) 

so we have the following transformation 

Q> = (2.8) 

Similarly, the momenta pi are Lagrange multipliers, so that their transformations are also 
not fixed a priori. We fix their transformations by requiring that the modified Lagrangians 
are the same: 

L*(q\ Q\ Q l ) = L*{cf>\ ¥, ¥) = L{(f>\ ¥) + - $ l ). (2.9) 

Then we have 

Pi(.q i -Q i )=nS i -&)- (2-10) 
Using the first equation of (2.7) and (2.8) in (2.10), we have 

* = %Vr (2.11) 

The Hamiltonian, obtained by the Legendre transformation of the right hand side of (2.9), 
is expressed as 

H* = mft + — L* — 7r, t ¥ + U, t ¥ - L. (2.12) 
Since L = L, the change of the Hamiltonian is 

AH* = H*-H* = IU& - PiQ\ (2.13) 
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Now 

n.^ = ^ = p.^ (214) 

and 

B 2 f i Ftp . 

*' = ad?i tQi + W Q '- (215) 

Using (2.15) and (2.16) in (2.13), we have 

AH " = wh QtQ ' n - (216) 

Thus the Hamiltonian changes under the transformation of the original generalized coordi- 
nates and reasonable transformation of other canonical variables. If the Hamiltonian repre- 
sents the energy of the system, which is often the case, this is unreasonable physically. The 
difference depends on the variables characteristic to higher derivative theory. We note that 
if the variables transform to make the Hamiltonian invariant, the Lagrangian is changed. In 
order to resolve this problem, we propose a new canonical formalism which can be interpreted 
to be the generalization of the Ostrogradski's method. 



3 A canonical formalism in terms of Lie derivatives 

The idea is simple. As noted above, the problem of Ostrogradski's method was resolved by 
the method of BL where the choice of the generalized coordinates corresponding to the time 
derivatives of the original generalized coordinates is extended to almost arbitrary functions 
of these variables. In the application to gravity, the extrinsic curvature is chosen as the 
new generalized coordinates. It is the Lie derivatives of the metric which is the generalized 
coordinates. Lie derivatives reduce to the time derivatives for flat space. So this choice 
can be thought of the simple generalization of the Ostrogradski's choice. Thus we move 
to the canonical formalism by the Ostrogradski transformation contrary to the Legendre 
transformation with the modified Lagrangian as in the method of BL. 

Ostrogdadski transformation is carried out as follows. Let us consider a system described 
by a Lagrangian containing the n-th time derivatives of the generalized coordinates q % (i = 
1 V). 

L = L(q\q\...,q^) (3.1) 
Take the following variation of the action 

5S = / L(q l + 5q\ . . . , q l{n) + 5q l{n) )dt - / L(q\ . . . , q l(n) )dt (3.2) 
Jti+Sti h x 

where 5q l are decomposed as the sum of the variations of the function q l (t), which we denote 
by 5*q\ 5*q l = (q+Sq)(t)—q(t), and those due to the change of the time , q(t+8t)—q(t) = q l St: 

5q' = 5*q l + q l bt (3.3) 

where the second terms are assumed to contribute only near the end points. Then the above 
variation is written, retaining only the first order terms in small quantities, as 

rti rt2+St 2 rti 

5S= L(q\...,q l ^)dt+ L(q 1 ,...,q t{n) )dt + 5* Ldt. (3.4) 

Jti+Sti Jt 2 Jti 
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Using the approximation 

ft k +St k 

the first two terms in (3.4) become the following 



L(q\ . . . , q^)dt = L (q\t k ), q l{n \t k )) St k , (k = 1, 2) 



[L5t\H = ^i. 



Rearranging the sum in the third term of (3.4), we have 
where 

^ (_l)r-«-l£>r- a -l£ 



TV 

<^ = E 

i=i 



5S 2 = 8* J* Ldt = [SF}\\ + <T;S 2 

<9L 



and 



s=0 kr=s+l 



TV ft2 » 



dL 



>( j*o<W 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



where represents the time derivatine, i.e. D = —. 5*S2 vanishes if we require the 

at 

variational principle. The new generalized coordinates are taken as 



ql = D s q i 



(3.10) 



and the momenta canonically conjugate to these cordinetes are defined by the coefficients of 
the variations of these coordinates in SF 



pl= E 

r=s+l 



^r—s—lj-jr—s—1 ^ 



dL 

dD r q i 



(3.11) 



The Hamiltonian is defined as (— l)x the coefficient of St in LSt + SF after using (3.3) 

(3.12) 



TV n-1 
i=l s=0 



Note that for s — n — 1, (3.11) has a simple expression 

dL 



V 



n-1 



dq l n -i 



(3.13) 



Thus the Ostrogradski transformation is a generalization of the Legendre transformation. It 
is noted that the highest order derivatives need not be the same for all i. 

In this work we investigate the higher-curvature gravity (HCG) of f(R) type in which the 
Lagrangian density is given by a function of the scalar curvature 



where g = detg^ v . The action is given as usual: 

S = J' 2 Ldt, L = J Cd 3 x 



(3.14) 



(3.15) 
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We adopt the ADM variables iV(x, t), iV l (x, t) and /iy(x, t) with respect to a hypersurface 
t = constant(E t ) as the generalized coordinates. In terms of these variables, the scalar 
curvature is expressed as follows 



R = 2N~ 1 h ij d K ij - U - 2N~ 1 AN - 2N~\N k d k K + 2N i > j K ij ). 



Here is the extrinsic curvature of E t 



(3.16) 



2N 



(d h i:j - N rJ - Nj.i) , 



and 



U 



3KijK ij - K 2 



R, 



(3.17) 



(3.18) 



where K = h 1 ^ and i? is the scalar curvature of S t . The determinant g is expressed as 
-iV 2 /i with h = det ft,. 



vi- 



lli this case <5F takes the following form 
dC _ „ ( dC \ 



d 6 x 



d(d h i:j 
dC 



-5*hij - d 



d(d N) 



5*N + 



d(dlhij) 
dC 



5*hij + 



dC 



d(d d k N* 



5*d k N l + 



d{dod k hij 
dL 



-5*d k hij + 



dC 



d{d 2 hij) 



5*d h. 



d(d N 



From this expression, it appears that N and N l , which represent the choice of the coordinate 
system, have the momenta canonically conjugate to them. In that case, they are allowed to 
be only the solutions of dynamical equations and the general covariance is broken. However, 
the time derivatives, d^h^, d N, dnN 1 and dnd k N l are involved only through d Kij. It is also 
noted that time derivatives d^hij and dnd k hij are involved through d^K^, and d k Kij. 
Thus 5F reduces to the integration of 



dC 



dC 



d(d h i:j ) 



dC 



dC 



\d{dlhij 
dK kl 



-d k 



dC 



d(d K kl ) dh 



+ d v 



V) 



d(d G d k hij) J 

' dC 8K kl ' 
d(d K kl ) d(d m hij) 



(3.19) 



8*K 



where we use the relations, e.g. 



d(d K, 



kij 



dK, 



ki 



d(d N) 



dN 



From this expression, if we adopt as the new generalized coordinate, instead of dohij, 
time derivatives of N and N l are absorbed in that of and the restriction on N and N l 
disappears. So is taken as the new generalized coordinate as in the method of BL and is 
denoted as Qij. The momenta canonically conjugate to hij and Qij, p 1 ^ and 11^ respectively, 
are taken as the coefficients of their variations in (3.19), according to Ostrogradski: 



:>:i 



V 



and 



dC 



d(d h 



dC 



d(d 2 h ij ) / 



-d k 



dC 



dC dK h 



kl 



d(d d k hij) 



d(d K k i) dhij 



f<9„ 



dC 



dK, 



ki 



d(d K k i) d(d m hij)J 
(3.20a) 



n 



dC 

d{doQij) ' 



(3.20b) 
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Hamiltonian H is defined as the coefficient of (—1) x St in L5t + 5F after using (3.3). Then 
Hamiltonian density is given by 



H = p ij d h tj + \Y ! <W,, - C. 



(3.21) 



It is noted that Kij is half the Lie derivative of along the normal to the hypersurface 
Yj t . The righthand sides of (3.20) are rewritten explicitely in terms of the Lie derivatives as 
follows 

4i 



p°-> = -y/h [f{R)Q ij + h ij f"{R)C n R + N' 2 d k NN k f'{R)h ij 

U ij = 2^hf\R)h ij 
where the scalar curvature is expressed as 

R = 2ti j C n Q ij + Q 2 - ^QijQ ij + R- 2AQnN). 

It is seen from (3.22) that IT 7 has only the trace part and is written 



n 
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\lih ij with n = 2dVhf'(R) 



d 



(3.22) 



(3.23) 



(3.24) 



where d is the dimension of space. Converting this relation, the scalar curvature is expressed 
as 

R = f'~\]l/2d\ / h) = ^(U/2dVh). (3.25) 

4 Invariance of the Hamiltonian 

In this section, we demonstrate that the Hamiltonian defined in the previous section is 
invariant under the transformation of the generalized coordinates induced by (1) the general 
coordinate transformation on the hypersurface S t and (2) the transformation of the metric 
hij to its function Gij(hki) such as the transformation of the scale factor a(t) to its logarithm 
4>{t) = lna(t). 



4.1 Invariance under the coordinate transformation on Y> t 

Let us consider the following general coordinate transformation on S t 



x 

X 1 



x° = x° 



x l 



f l (x l ,x 2 , X 3 ). 



(3.26) 



Under this transformation, transforms as a second rank tensor, iV a scalar and N l a 
vector. On the other hand the momenta transform as contravariant tensor densities 



d(x) dx l dx j kl 
d(x) dx k dx l ^ 

d(x) dx* dx j u ki 
d(x) dx k dx l 



(3.27) 
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where a overbar represents the transformed quantity and d{x)/d{x) is the Jacobian of the 
transformation. The transformed Hamiltonian density is defined as in (3.21) 

U = p'dohij + rFdoQy - C. (3.28) 

The time derivatives and the Lie derivatives do not effect the transformation properties so 
that Qij, dohij and d^Qij transform as second rank tensors. Since the Lagrangian density is 
the scalar density, so is also the Hamiltonian density 

o(x) 

Therefore we have the desired result 

H = t Hdtd 3 x = I'" ^lndt^-d 3 x = H. (3.29) 
Jti Jti o(x) d(x) 

This result is expected since the coordinate system we begin is not specified. Thus what we 
have shown is essentially the consistency of our method. 



4.2 Invariance under the transformation to the function of the 
metric 



Next consider the transformation of the generalized coordinates 

hij — * <f>ij = Fij(hki) or hij = Gij{4>ki). 



(3.30) 



Under this transformation, three dimensional space is unchanged so that the new generalized 
coordinate Qij = C n hij/2 is unchanged and so is the momentum canonically conjugate to it 
IP- 7 = dC/d(doQij). These are expressed in terms of the transformed quantities as 



, IF = 2Vhf^)G^{(t> kl ). 
Then using (3.22), (3.31) and the relation 

dGi, 



(3.31) 



Shij = 

.90 



(3.32) 



the first term of SF, (3.19), takes the following form 



-Vh{f'(iP)Q^ + Gtif'WCrd + N- 2 d k NN k f'(R)h"} ~Q^f<t>ki + 2y/hf'W)G*5*Q ij . 

(3.33) 

The momentum canonically conjugate to 0^ is defined as the coefficient of 5*0^: 



P? = ~Vh{f'WQ kl + G kl f'^)C n ^ + N- 2 d k NN k f{R)W} 



dG 



ki 



'13 



(3.34) 
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From (3.22) and (3.33), we have the following relation 




(3.35) 



Transformed Hamiltonian density Ti^ is defined as 



ft = p;^.. + rFQ 



c. 



(3.36) 



The change of the Hamiltonian density is given by 



An = n^-n = v 13 ^, 




(3.37) 



Using (3.32) and (3.34), it is shown that AH. vanishes. This result should be necessary if 
the Hamiltonian has something to do with the energy. 

Application to the FRW spacetime, in which case the transformation is from the scale 
factor a{t) to its function G(a), e.g., In a(t), is straightforward. 

5 Summary and discussions 

We proposed a canonical formalism of HCG by combining those of Ostrodradski and BL 
using the Lie derivatives instead of the time derivatives. Invariance of the Hamiltonian 
is shown under the transformation of generelized coordinates preserving the hypersurface 
£j which is lacking in BL. In fact, the transformation properties in HCG have not been 
addressed. 

The result is important in HCG where the conformal transformation to Einstein frame 
is an often used technique [18]. The problem with this transformation is which of the metric 
is physical, i.e. the observed one. Many criteria have been proposed to select the physical 
one, however the problem has not been settled[19]. However the conformal transformation 
depends on the curvature, so does on the momenta from the viewpoint of canonical formalism. 
It is not obvious whether the canonical equations of motion in both frame is equivalent. Also 
it would be necessary to show that the Poisson brackets defined in each frame are consistent. 
If the answers to these problems are negative, equivalence of both frames is broken and the 
canonical quantization leads to different quautum theories. The results of these problems 
will be reported separately. 

Finally we comment on quantum gravity. It is well known that there are two versions. 
One is the quantum theory of gravitons based on the duality of gravitational wave and 
graviton and has the same footing as the ordinary quantum theory of matter and radiation. 
This version, however, requires a background spacetime. If the quantum theory is the fun- 
damental theory in physics, the background spacetime should also be determined quantum 
theoretically. It is the subject of quantum cosmology which is the second version of quantum 
gravity. In this version no guiding principle such as duality of particles and waves has been 
known. The usually adopted procedure is the canonical quantization. For this, the classical 
theory is a prerequisite. Such a theory may not be the Einstein gravity but a kind of HCG. 
If this were the case, investigations of the early stage of the univerese, or string theory might 
give some clues. The strings feel tidal force, which necessarily leads to HCG. Observational 
informations are generally more important. These informations may be brought about from 
future cosmological observations. 
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